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Abstract. A widely applied method of supported lead estimation in sediments using 
gamma spectrometric *?°Ra determination via *”*Rn short lived daughter products 
relies on radioactive equilibrium between Ra and *”*Rn being established after 
sealing the samples. Advantages and disadvantages of methods of ?7°Ra estimation in 
sediments, using either *?°Ra daughter products or direct estimation by the 186.2 
keV gamma-emissions are discussed. An equilibrium experiment was performed 
using test samples and in one case radioactive equilibrium was not reached. On 
theoretical sediment profiles it is shown how systematic errors in supported *!°Pb 
estimation can lead to wrong interpretations of ?'°Pb,, profiles and therefore affect 
210Pb derived chronologies. 


1. INTRODUCTION 


*1°Pbh sediment chronology is a frequently used tool for dating recent lacustrine 
and marine sediment cores or peat bog accumulations [1]. Its range is 
approximately the last 100-150 years, which is the period of the most significant 
anthropogenic effects to the environment. 7!°Pb (half-life 22.3 years) in 
sediments consists of two components: 

supported *'°Pb, which is present due to autigenic material of the 

sediment and 

unsupported (or excess) *!°Pb, which originates from the atmospheric 

deposition. 
Excess *!°Pb is then determined by subtracting supported activity from the total 
activity and used for estimation of accumulation rates and age models. 
Measurement techniques of total signal of 7!°Pb include either alpha 
spectroscopy via 7!°Po with assumption of radioactive equilibrium between ?!°Po 
and ?7!°Pb, or gamma spectroscopy. Both approaches have certain advantages 
and disadvantages [2]. In the authors’ laboratory, gamma spectroscopy is used. 
That allows simultaneous estimation of other natural and artificial gamma 
emitters within a single spectrum and the analysis is non-destructive. 


2. SUPPORTED 2'°PB ESTIMATION 


For supported ?!°Pb, which is usually assumed to be in radioactive equilibrium 
with its parent nuclide 77°Ra (half-life 1600 yr), analytical options were 
suggested by Moser [3]. The simplest method, sometimes used in combination 
with alpha spectroscopic determination of 7'°Pb, is to use the lowest constant 
activity of a sediment profile (samples with an age greater than 100 years). This 
method assumes a constant activity of **°Ra over the entire sediment profile, 
which is not always the case. It is also possible to measure Ra by alpha counting 
or liquid scintillation after radiochemical separation, or gamma spectroscopy [3, 
4]. 


2.1. Using ??°Ra daughter products 


Gamma spectroscopy of daughter products of **°Ra, namely short lived *‘*Pb and 
214Bi (further on named “daughter products”), is a widely used method for ?®Ra 
determination. Here we have to deal mainly with difficulties related to *?*Rn 
diffusion and recoil emanations from the sample, leading to systematic 
underestimation of Ra concentration [5]. In the study of Stoulos et al. [6] Rn 
gas leaking from the sample container, its accumulation in the void upper part 
of the container and attachment of Rn decay products onto its inner surface, 
resulting in another geometry, and self-absorption from that of the calibration 
sample, was resolved by adding active charcoal to the sample. However, this 
would change properties of the samples, which contradicts one major advantage 
of gamma spectroscopy for sediment chronology being a non-destructive 
method. An optional solution is using Rn proof steel containers [6], which would 
lead to decline in sensitivity in detection of low energy gamma radiation of *!°Pb. 
A common practice in many gamma spectroscopic laboratories is sealing 
samples using gas-tight foil and waiting for at least 21 days to allow ingrowth of 
daughter products. 


2.2. Direct estimation 


226Ra can also be directly estimated by its gamma line at 186.2 keV, but 
interference with the most prominent line of **°U 185.7 keV comes into account. 
A method for correction of the *°U and ??°Ra interference using the 63.3 keV 
line of **Th as a proxy for *°U was proposed by Dowdall et al. [7]. In this case 
constant natural ratio to ?U is assumed. When using the 92.6 keV doublet of 
*34Th, precautions must be taken due to its potential interference with Ka line 
from *°*Th. Contribution of *°U in the 186 keV peak can also be determined 
using another less intensive *°°U line at 143.8 keV [8], if detectable. 


2.3. Possible 22°Ra underestimation 


Motivation of this study was the participation of the laboratory in a laboratory 
intercomparison. It resulted in recognition of the possibility of underestimating 
226Ra values. It indicated that the existing laboratory practice of estimation of 
26Ra via daughter products using plastic sample dishes sealed in assumed gas- 
proof metallic foil, was possibly a weak point. As a countermeasure, the 
laboratory adopted new thicker metallized plastic foil-bags for sealing samples 
and performed an experiment using soldered steel containers. 


3. MATERIALS AND METHODS 
3.1. Test samples 


The test sample, denominated RV VII-2006, was a soil material used for the 7th 
Intercomparison of environmental samples organized by the German Federal 
office for Radiation Protection (BfS) [5]. It was taken from an area in federal 
state of Thuringia with naturally increased **°Ra levels. The test material had 
been grounded into a fine powder, so that 98.8 % of the fraction had a particle 
size less than 315 um. A cylindrical plastic dish with diameter of 70 mm and 
height of 20 mm was filled with the test sample. It was then sealed in a 
metallized plastic foil, considered to be gas-tight. Another experiment was 
performed using a steel container with diameter of 101.5 mm and height of 31.5 
mm in which the test sample was soldered. 


3.2. Gamma spectroscopy 


The samples were measured by low-level low-background gamma spectroscopy 
using a coaxial HPGe detector Canberra Industries (51% rel. efficiency) housed 
in a 10 cm Pb shielding with Cu, Cd and plastic lining. The spectra were 
analyzed by Canberra Genie 2000 software. Efficiency calibration was 
performed using the Monte Carlo based LabSOCS calibration tool [9], which 
takes into account not only sample to detector geometry, but also sample 
density and composition, as well as measurement container properties. 


4. Experimental, results and discussion 
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Figure 1. Establishing of radioactive equilibrium in test sample RV VII-2006 within the 
time interval of approximately 8 weeks for (a) a plastic dish sealed in a metallized 
plastic foil and (b) a soldered steel container. In the first case, the *!*Pb values reached 
only 88% of that of parent 7*°Ra, which was determined by 2 different methods (1) and 
(2) listed in chapter 4.1. In the steel container the radioactive equilibrium was reached. 
The dashed line indicates the reference value. The error bars for all data are 1 standard 
deviation (they include counting statistics and calibration of the measuring device). 


4.1. Equilibrium test 


Gamma spectra of the BfS test samples sealed in 2 upper mentioned ways were 
recorded repeatedly during approximately 8 weeks (16 ???Rn half-lives) following 
their sealing in order to observe the process of radioactive equilibrium 
establishing. The daughter product **Pb was measured via its 242.0, 295.2 keV 
and 351.9 keV gamma lines. ?*°Ra activities were determined using the gamma 
line 186.2 keV after subtraction of the *°U contribution, for which 3 methods 
were applied: 
(1) using the *°°U line at 143.8 keV (after interference correction from 
contributing °?” Ra, determined by the 144.2, 154.2, 269.5, 323.9 and 
338.3 keV lines), 
(2) using 63.3 keV and 92.6 keV line of Th and 1001.0 keV line of 
*34mPa for estimation of *°8U with activities derived from the natural 
238U/235U activity ratio of 21.7. 


These both methods applied on the metallized plastic foil sealed sample lead to 
very consistent results with no statistical difference (two sample t-test 
(34)=0.46, p=0.65) (Fig. 1a). 

(3) Finally, a combination of methods (1) and (2) was used, applying all 4 
mentioned lines for ?U contribution estimation, which resulted in 
reduction of ?®Ra uncertainty by 43%. 

While the mean value of Ra measured by all 3 methods was close to the 
reported reference value, the equilibrium value of ?'*Pb in the metallized plastic 
foil sealed sample reached only 88% that of ”°Ra. The mechanism of this 
phenomenon, which has not been explained within this study yet, is speculated 
to be Rn diffusion through the foil. The very fine material used for the test is 
prone to show a high emanation rate. Efficiency problems with daughter 
product estimation can be excluded, because both daughter products *!*Pb and 
214Bi gave consistent results when estimated by lines of wide energy span 
(between 242 and 1764 keV). In the parallel experiment using soldered steel 
container the radioactive equilibrium was reached (Fig. 1b). 


4.2. Effects of supported 7!°Pb underestimation on excess 7!°Pb 
sediment depth profiles 
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Figure 2. Theoretical ?Pb sediment depth profiles described in chapter 4.2.: (a) a 
simple profile (Profile A) and (b) a profile with a random scatter (Profile B). For both 
profiles, in the left figures, °Pbi and both correct and systematically underestimated 
2°Pbsup are depicted. In the central figures, the correct and overestimated 7!°Pb,; 


(calculated as 7!°Pbyot-*'°Pbsup), together with simple 2 parameter exponential fits, using 
least square algorithm, are plotted. In the right figures, residuals to the fits are shown. 


On 2 theoretical models it is illustrated how systematic underestimation of 
21°Pb.uy» Can affect the resulting sedimentation rate estimation in *'°Pb,, sediment 
chronology. The simple example (Profile A) (Fig. 2a) presents a theoretical non- 
compacted and non-mixed sediment profile with a constant sedimentation rate 
of 0.25 cm-y’. The sampling interval is 1 cm and the depth of the profile is 30 
cm. 7!°Pbsup being constant 25% of the total ?'°Pb initial activity (top of the 
profile) is subjected to underestimation by 12%. The value of 25% for ratio 
between supported lead and total initial 7‘°Pb was selected as a “typical value”, 
although based on our experience, it can vary strongly according to local 
environment, often between 5-50%. Underestimation of supported lead in this 
example leads to a constant shift of 7!°Pb,, to greater values. 

Another example (Profile B) shows the upper described model modified by 
adding a random scatter (using normally distributed random number generator) 
to ?Pbio and 7!°Pb.., data and error bars added (Fig. 2b). 

On a real dataset (with scatter caused by counting statistics and other natural 
causes) excess 7!°Pb overestimation can remain unobserved, and could also be 
attributed to not reaching the “dating horizon”. It can possibly be detected by 
careful examination and finding regular patterns in the residuals plot (after 
fitting) (Fig. 2a), but this effect is not obvious in naturally scattered datasets 
(Fig. 2b). The described effect will be strongest in environments with higher 
proportion of supported 7!°Pb. Systematic error in *'!°Pbsu» estimation can 
therefore have consequences for *!°Pb,, derived chronologies leading to 
systematic underestimation of ages. 


4.3. Possible consequences of supported ?!°Pb underestimation for 
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Figure 3. Chronologies derived from the 2 theoretical profiles: (a) Profile A and (b) 
scattered Profile B. For both cases chronologies were derived by both CIC and CRS 


models using both correct and overestimated ?'°Pb,, profiles. 


Two most commonly used chronology models were applied to the upper 
described 2 theoretical profiles A and B: Constant initial concentration (CIC) 
model [10] and Constant rate of supply (CRS) model [11]. 

In the first case (CIC model), 7!°Pb,, profiles were fitted by a simple 2 parameter 
exponential function using least square algorithm. Based on resulting 
parameters, sedimentation rates were calculated. The resulting sedimentation 
rates are 0.284 + 0.004 cm-y" for Profile A and 0.280 + 0.026 cm-y' for Profile 
B. Sedimentation rates were thus overestimated by over 13%. 

For comparison, 7'°Pb,, profiles were used for calculation of ages according to 
CRS chronology model, which in case of non-scattered Profile A led to even 
more serious age underestimation than CIC model. Graphical depth/age 
relationship for both CIC and CRS models for Profiles A and B are presented in 
Fig. 3. 


4.4. Discussion 


A method for testing and correcting possible underestimation of 7!°Pb.u» in 
sediment profiles measured by gamma-spectroscopy is currently being 
developed. In case of supported lead underestimation caused by partial Rn 
escape from the measurement containers, corrections can be employed. 


5. CONCLUSIONS 


It was shown that although samples for *”°Ra analysis are carefully sealed in 
presumably gas tight foils, especially in very fine grained samples the 
radioactive equilibrium required for **°Ra estimation via daughter products 
might not be reached. 

This can theoretically have consequences also for supported 7!°Pb estimation via 
the Ra daughter products. A systematic underestimation of 7!°Pbsup can result 
in systematic underestimation in 7'°Pb,, derived ages, as it was presented on 
theoretical examples. 
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